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Abstract 

Solutions of an implicit ODE form a web. Already for cubic ODEs the 3-web of solutions 
has a nontrivial local invariant, namely the curvature form. Thus any local classification 
of implicit ODEs necessarily has functional moduli if no restriction on the class of ODEs is 
imposed. Here the most symmetric case of hexagonal 3-web of solutions is discussed, i.e. the 
curvature is supposed to vanish identically. A finite list of normal forms is established under 
some natural regularity assumptions. Geometrical meaning of these assumptions is that the 
surface, defined by ODE in the space of 1-jets, is smooth as well as the criminant, which is 
the critical set of this surface projection to the plane. 

Key words: implicit ODE, hexagonal 3-web, equivariant diffeomorphism. 
AMS Subject classification: 37C15 (primary), 53A60, 37C80 (secondary). 

1 Introduction 

Consider an implicit ordinary differential equation 

F(x,y,p)=0 (1) 

with a smooth or real analytic F. This ODE defines a surface M: 

M := {(x,y,p) Gl 2 x P 1 (M) : F(x,y,p) = 0}, (2) 

where (x,y,p) are coordinates in the jet space J 1 (M, R) with p = Generically the condition 
grad(F)\ F ^ x y p ^ =0 / holds true for any point m = (x,y,p) £ M, i.e. M is smooth. If the 
projection ir : M — » M?, (x,y,p) i— > (x,y) is a local diffeomorphism at a point m £ M then this 
point is called regular. In some neighborhood of the projection of a regular point 7r(m) equation 
(HJ) can be solved for p thus defining an explicit ODE. 

If the projection ir is not a local diffeomorphism at m, then the point m is called a singular 
point of implicit ODE ([1]). The set of all singular points is called the criminant of equation (P) 
or the apparent contour of the surface M and will be denoted by C: 

C := {(x,y,p) El 2 x F 1 ^) : F(x,y,p) = FJx,y,z) = 0}, (3) 
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where the low subscript denotes a partial derivative: F p = 

Studying of generic singular points of implicit ODEs was initiated by Thorn in [23] . Due to 
Whitney's Theorem such points are folds and cusps of the projection ir. Local normal forms for 
generic singularities were conjectured by Dara in [9] and for a generic fold point were established 
by Davydov in [10]. The classification list for a generic fold point of the projection ir is exhausted 
by a well folded saddle point, a well folded node point, a well folded focus point and the regular 
singular point, where contact plane is transverse to the criminant. Cusp points were studied by 
Dara [9], Bruce [5] and Hayakawa, Ishikawa, Izumiya, Yamaguchi [18J. 
Usually the following regularity condition is imposed at each point of the criminant: 

rank((x,y,p) ^ (F,F P )) = 2. (4) 

This regularity condition implies that the criminant is a smooth curve. At each point m outside 
the criminant C the contact plane dy — pdx = cuts the tangent plane T m M along a line thus 
giving a direction field, which takes the form: 

T=[F p :pF p :-(F x +pF y )], (5) 

in the coordinates (x,y,p). This direction field is called the characteristic field of M. The 
projection 7r(7) of an integral curve 7 C M of the characteristic field r is called a solution of 
ODE dl]). If F pp 7^ at a point m £ C on the criminant then ([1]) reduces locally to quadratic 
in p (or binary) ODE. Such equations were the subject of intensive study. See, for example, 

Suppose equation ([1]) has a triple root po at (xo, yo) then the equation F = can be written 
locally as a cubic equation 

p 3 + a(x, y)p 2 + b(x, y)p + c(x, y) = 0, (6) 

as follows from the Division Theorem. Thus if in a domain U C M? outside the discriminant 
curve A := ir(C) this cubic equation has 3 real roots Pi,P2,P3, we have 3- web formed by solutions 
of ([1]). A generic 3- web has a nontrivial invariant. In differential-geometric conext this invariant 
is the curvature form of the web. Therefore any general local classification of cubic implicit 
ODEs ([6]) necessarily has functional moduli (cf. TOj). Moreover, this invariant is topological 
in nature hence even topological classification will have functional moduli if no restriction is 
imposed on the class of ODE. (See also [20] and [21], where web structure was used for studying 
geometric properties of differential equations.) 

In this paper we consider cubic ODEs ([6]) with a hexagonal web of solutions. Equations of 
this type describe, for example, webs of characteristics on solutions of integrable systems of three 
PDEs of hydrodynamic type (see [12], [13]). Another example is WDVV associativity equation 
(see Example [2] below) . 

Definition 1 Let U C I 2 5e the open set, where (0|) has 3 real roots pi,P2,P3 o- n d suppose 
U 7^ 0. We say that ([6]) has a hexagonal 3-web of solutions if for the projection ir{m) of each 
regular point m € M with ir{m) 6 U there is a local diffeomorphism at ir(m) mapping the 
solutions of (0) to three families of parallel lines. 

The case of hexagonal web of solutions is also the most symmetric, i.e. the Lie symmetry 
pseudogroup of ([1]) at a regular point has the largest possible dimension 3. The list of normal 
forms turnes out to be finite provided regularity condition Q is satisfied. These forms are given 
by the following examples. 
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Example 1 The classical Graf and Sauer theorem [TB] claims that a 3-web of straight lines is 
hexagonal iff the web lines are tangents to an algebraic curve of class 3, i.e. the dual curve 
is cubic. This implies immediately that the following cubic Clairaut equation has a hexagonal 
3-web of solutions: 

p 3 + px — y = 0. (7) 

The solutions are the lines p = const enveloping a semicubic parabola. (See Fig. [T]) Note that 
the contact plane is tangent to M along the criminant, i.e. the criminant is a Legendrian curve. 

Example 2 Consider associativity equation 

Uxxx — ^xyy ^XXyUyyy, 

describing 3-dimensional Frobenius manifolds (see [UJ). Each of its solutions u(x,y) defines a 
characteristic web in the plane, which is hexagonal as was shown by Ferapontov [14]. Charac- 
teristics are integral curves of the vector field 

d x - X(x,y)d y , 

where A satisfy the characteristic equation 

A ~\~ HyyyX ^UxyyX ~\~ l^xxy — 0* 

For the solution u = — h ^ the characteristic equation becomes 

p 3 + 2xp + y = (8) 

after the substitution x — ► —x, y — > —y, A — > —p. The criminant of this ODE is not Legendrian 
and the solutions have ordinary cusps on the discriminant (see Fig. [1]). The discriminant is also 
a solution. In the analytic setting the above two normal forms were conjectured by Nakai in 




Figure 1: Solutions of p 3 + px — y = (left) and p 3 + 2xp + y = (right) with horizontal y-axis. 

We find also local normal forms at points, where the projection n has a fold, i.e. the cubic ODE 
factors out to a quadratic and a linear terms. 
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Example 3 Suppose the criminant of ODE is Legendrian then this ODE is locally equivalent 
to 

P 2 = V- (9) 

Solutions of this ODE together with the lines dx = form a hexagonal 3- web (see Fig. [1]). In 
fact, both the lines dy = or the parabolas 2dy — xdx = also supplement the 2- web of solutions 
of © to a hexagonal 3-web, but the surfaces of the corresponding cubic equations p(p 2 — y) = 
and (2p — x)(p 2 — y) = are not smooth at m = (0,0,0). If we agree to consider a quadratic 
equation as a cubic with one root at infinity, than equation ([9]) is the third normal form in our 
list. The following coordinate change y = y + x = x straightens the solutions, transforming 
ODE ([9]) to a quadratic Clairaut equation 

p 2 + px — y = 0. 

As the the lines dx = are preserved this example is also a special case of the Graf and Sauer 
Theorem. 

Example 4 Suppose the criminant of ODE is not Legendrian then this ODE is locally equiva- 
lent to 

p 2 = x. (10) 

Solutions of this ODE together with the lines dx = form a hexagonal web (See Fig. [2]). The 
lines dy = also complete the 2-web of solutions of (jlOp to a hexagonal 3-web, but again the 
surface M of the corresponding cubic equation p{p 2 — x) = is not smooth at m = (0, 0, 0). 




Figure 2: Solutions of p 2 = y and the lines dx = (left). Solutions of p 2 = x and the lines 
dx = with horizontal y-axis (right). 



Example 5 For completeness let us mention the case of a regular point of an implicit cubic 
ODE. If its 3-web of solutions is hexagonal it can be mapped to the web of 3 families of parallel 
lines dx = 0, dy = and dx + dy = 0. This gives the following "cubic" ODE 

p(p +1) = 0. (11) 

Now we can formulate our classification theorem. 
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Theorem 1 Suppose functions a, b, c are real analytic and the following conditions hold for an 
implicit cubic ODE 



F(x, y,p) := p 3 + a(x, y)p 2 + b(x, y)p + c(x, y) = 

at a point m = (x ,y ,po) G M := {(x,y,p) £l 2 x P 1 (R) : F(x,y,p) = 0}: 

1 ) this equation has a hexagonal 3-web of solutions, 

2) dF\ m / 0, 

3) rank((x, y,p) i— > (.F, i^,))| m = 2 if m lies on the criminant C. 

Then this ODE is equivalent to one of the following five forms with respect to some local real 
analytic isomorphism: 

i) p 3 + 2xp + y = 0, if po is a triple root and the criminant is transverse to the 

contact plane field in some punctured neighborhood of m, 

ii) p 3 + px — y = 0, if Po is a triple root and the criminant is Legendrian, 

Hi) p 1 = y, if Po is a double root and the criminant is Legendrian, (12) 

iv) p 2 = x, if Po is a double root and the criminant is transverse 

to the contact plane at m, 

v) p(p + 1) = 0, if the roots are pairwise distinct at n(m) = (xo,yo)- 

If the functions a,b,c are smooth and conditions 1),2),3) are satisfied then there is a diffemor- 
phism of a neighborhood of the point (xQ,yo) onto a neighborhood of the point (0,0) reducing 
the above cubic ODE either to one of the four equations (T2Hi)- (TEv) or to an equations that 
coincides with within the domain, where hl2i) has three real roots. 

The main difficulty in proving the above classification theorem brings the case of irreducible cu- 
bic ODE. The idea is to lift its 3-web of solutions to M and then to the plane E : p\ +p2 +P3 = 
in the space of roots of the cubic equation p 3 + A(x,y)p + B(x,y) = 0. (Note that the gen- 
eral case reduces to this cubic.) Then this 3-web at the plane E has D3-symmetry permuting 
the roots. Using the regularity condition we construct a D3-equivariant diffeomorphism "up- 
stairs" , matching the web to that of a corresponding normal form. Due to the D3-symmetry the 
constructed diffeomorphism is lowerable to some diffeomorphism "downstairs", i.e. to a point 
transformation in the plane of solutions. Most of the claims and the proofs below are given for 
the smooth case and for some neighborhood of the projection ir(m) of m G M, if it is not stated 
explicitly. In section [J] we discuss how to get rid of the annoying stipulation in Theorem [1] for 
the smooth case (|12i ) by replacing Definition [T] with a less geometric one. 



2 Normal forms for a fold point 

In this section we establish normal forms for the case, when the projection tt has a fold point 
at m. If cubic equation ([6]) has two coinciding roots then the the third root defines a regular 
point of the projection tt and the equation factors to a quadratic equation and a linear one. 
Regularity condition (Jl|) for a double root po implies immediately that the projection tt has a 
fold point at m. First we find a normal form for fold points and the symmetries of this normal 
form. Further we look for the linear in p (i.e. explicit) equation whose solutions complete the 
2-web of solutions of the quadratic normal form to a hexagonal 3-web. Finally, we bring these 
linear terms to some normal forms using the symmetries of the quadratic equation. We start 
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with a Legendrian criminant, then consider non-Legendrian criminant and finally show that the 
case of an isolated point of tangency of the criminant and the contact plane is excluded by the 
regularity conditions. 

2.1 The case of Legendrian criminant 

Proposition 1 Consider implicit ODE ([!]) with a smooth Legendrian criminant and a smooth 
surface Then characteristic field |5|) can be smoothly extended to the criminant. Moreover, 
the extended characteristic field is transverse to the criminant. 

Proof: Let m be a point on the criminant. A suitable contactomorphism ip maps M to M' with 
the following properties: 

a) the criminant of M is mapped to the line x = y = 0, 

b) <p(m) = (0,0,0), 

c) the tangent plane T m M is mapped to the plane y = 0. 

It suffices to prove the proposition for the transformed surface M' := {(x,y,p) £ M 2 x P 1 (M) : 
G(x, y,p) = 0} with the Legendrian line x = y = 0. Condition c) allows to rewrite G(x, y,p) = 
as 

y = g(x,p), 

while condition a) implies g(x,p) = xf(x,p). Now g p (0, 0) = g x (0,0) = implies /(0,0) = 
0, so that f(x,p) = pu(x,p) + xv(x,p) by the Hadamard lemma. As the line x = y = 
is Legendrian, the form dy — pdx must vanish on it: {d(xpu(x,p) + x 2 v(x,p)) — pdx}\ x=Q = 
{(pu(x,p)dx — pdx}\ x=Q = or n(0,p) = 1. Again by the Hadamard lemma one gets u(x,p) = 
1 + xw(x,p) and 

y = xp + x 2 h(x,p) 

with h(x,p) = pw(x,p) + v(x,p). Now the characteristic field r is defined by restriction of 
equation dy — pdx = to M', i.e. by {d{xp + x 2 h(x,p)) — pdx}\ M , = x{(l + xh p {x,p))dp + 
(2h(x,p) + xh x (x,p))dx} = 0. This implies that the characteristic field on M' in coordinates 
(x,p) is generated by the vector field (1 + xh p (x,p))d x — (2h(x,p) + xh x (x,p))d p which is clearly 
smooth and transverse to the line x = y = on M' . 

□ 

Theorem 2 Let (OP be an implicit ODE such that corresponding surface (dj) is smooth. Suppose 
its criminant C is a smooth Legendrian curve and the projection n : (x,y,p) — » (x,y) has a fold 
singularity at m G C. Then ([!]) is locally equivalent to 

p 2 = y (13) 

with respect to some diffeomorphism ip : U — > U , where U,U £ E? are neighborhoods of (0,0), 
ir(m) and </?(0,0) = 7r(m). 

Proof: Due to Proposition [1] there is a smooth direction field r on M. Its integral curves define 
a foliation J 7 of a neighborhood of m. As the projection tt has a fold on the criminant, the 
discriminant curve A = tt(C) G M 2 is smooth. Let us choose new coordinates on (U,ir(m)) 
such that the discriminant curve turns to the line y = 0. Then equation ([T]) is equivalent to 
y = g(x,p). The discriminant curve is a solution therefore g(x,p) = pf(x,p) by the Hadamard 
lemma. The criminant y = p = is Legendrian hence {d(pf(x,p)) —pdx} p= o = f(x, 0)dp = or 
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f(x,0) = 0. Now f(x,p) = ph(x,p) and y = p 2 h(x,p). As the projection ir has a fold at m £ M 
holds true 9p(p 2 /i(x,p))| a; Q = 2/i(0,0) / 0. Consequently characteristic field © on M with 

F(x,y,p) = p 2 h{x,p) — y is generated by the vector field 

(2ph(x,p) +p 2 h p (x,p))d x - (p 2 h x (x,p) - p)d p = p{(2h(x,p) + ph p (x,p))d x + (1 - ph x (x,p))d p }. 

Due to /i(0, 0) ^ it is transverse to the kernel of dir on M hence the projection of each integral 
curve crossing the criminant C is smooth and tangent to the discriminant curve. Locally the 
equation p 2 h(x,p) — y = can be rewritten as a quadratic equation 

p 2 + a(x, y)p + b(x, y) = 0. 

This easily follows from the Division Theorem since h(0, 0) 7^ 0. Moreover, as the criminant 
y = p = is a Legendrian curve holds true b(x, 0) = a(x, 0) = 0. Thus one gets a(x, y) = ya(x, y) 
and b(x,y) = y0(x,y) with /?(0,0) 7^ as M is smooth at m = (0,0,0). Consider (x,p) as local 
coordinates on M and define the map i : M — > M by 

i(x,p) = (x, —p — p 2 h(x,p)a(x,p 2 h(x,p))). (14) 

This map is the involution that permutes the roots pi,P2 of our quadratic ODE. Let the foliation 
T on M be defined locally by I(x,p) = const with grad(/)| a , =p=0 7^ 0, where / is a first integral 
of the characteristic field r. Then the functions I and J := are functionally independent as 
di{r) and r are transverse to each other. Since r is transverse also to the criminant the partial 
derivative I x \ x=p=0 does not vanish. Hence one can chose / so that I(x, 0) = x which implies 
J(x, 0) = x. Let us take thus normalized functions /, J as local coordinates on M. Note that 
the following relation holds true 

tt(I,J)=tt(J,I), (15) 

since i{I,J) = (J, I)- For the ODE p 2 = y the above defined objects are as follows: I(x,p) = 
x — 2p, i(x,p) = (x, —p), J(x,p) = x + 2p. Now define a diffeomorphism germ if) : M, — ► M, 
by ip(I,J) = (I, J)- By () 15f> there exists such a map germ 99 : M 2 ,0 -» R 2 ,0 that the following 
diagram commutes: 

M, ^- M, 

I vr i 7f (16) 

We claim that the map germ 99 is the searched for diffeomorphism. By construction it maps 
solutions to solutions. To show that ip is differentiable consider the map germ n o -0 : M, — > 
M 2 ,0. Applying Malgrange's Preparation Theorem to this map germ in coordinates (x,p) on 
M, and (x, y) on IR 2 , one gets 

x = X (x,p 2 ) +pX 1 (x,p 2 ), y = Y (x,p 2 ) +pY 1 (x,p 2 ) J 

where the functions Xq, X\, Yq, Y\ are smooth. Since ir o ifi(x,p) = it o ip(x, —p) = (x, y) holds 
true Xi(x,p 2 ) = Yi(x,p 2 ) = 0. With p 2 = y the above formulas take the form 

x = X (x,y), y = Y (x,y), 

and therefore define a smooth map germ ip. Applying the same considerations to the map germ 
tt o ip' 1 we see that ip has an inverse for y > 0. Therefore ip is invertible. □ 
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Remark. Note that equation (|13|) has a nontrivial symmetry group. For example, the scaling 
x — ► Xx, y —* X 2 y leaves it invariant. Therefore the diffeomorphism ip in Theorem [2] is not 
unique. 

Proposition 2 Equation [13\) has an infinite symmetry pseudogroup. Its transformations are 
given by 

x = F(x + 2^)+F(x-2^), y = h i F{x + 2^)-F{x-2^)} 2 . (17) 

Here F is a smooth function subjected to F'(u)\ u=Q ^ 0. Infinitesimal generators of this pseu- 
dogroup have the form 

{f(x + 2^y) + f(x - 2^y)}d x + Jy{f{x + 2Jy) - f(x - 2^)}d y , (18) 

where f is an arbitrary smooth function. 

Proof: Consider the action of symmetry group transformation on M in coordinates (u, v), where 

u = x — 2p, v = x + 2p. (19) 

To preserve the foliations T and i{T) determined by the direction fields r and di(r) it must have 
the following form u = 2F(u), v = 2G(v). This transformation is a symmetry if it commutes 
with i. Note that i permutes u and v and sends d p to — d p . Hence F = G. Now substitutions 
x = \(F(v) + F{u)) and y = p 2 = \{F(v) + F{u)) 2 gives (HH). Condition F'(u)\ u=0 / is 

equivalent to non-vanishing of the Jacobian: 

1 x=y=0 

Consider now infinitesimal symmetries of (|13j) . They are defined by operators £(x,y)d x + 
rj(x,y)d y . These operators must be liftable to M. On M the lifted vector field must be a 
symmetry of foliations T and i{T\ In coordinates (u,v) any infinitesimal symmetry X of 
foliations T and is easy to write down: 

X = f{u)d u + g(v)d v . (20) 

In coordinates (x,p) on M this operator X takes the form 



A = \{g(v) + f(u))d x + \{g(v) - f(u))d p . 



X is lowerable iif di{X) = X. Now lowerability condition amounts to f(u) + g(v) = f(v) + g{u) 
and g(v) — f(u) = f(v) — g(u). Thus the first equation gives f(u) — g(u) = c = const and 
the second implies c = 0. Substitution u = x — 2^/y, v = x + 2^/y, p = y/y into X gives 
(|18p for y > 0. The obtained transformations are correctly defined for y > and are smoothly 
(analytically) extendable for y > 0. The possibility to extend them for y < easily follows from 
Malgranges's Preparation Theorem: one considers the projection ir : M — > M. 2 , (x,p) t— > (x,p 2 ) 
in coordinates x,p on M and x, y on M 2 and observes that F(x + 2p) + F(x — 2p) is an even and 
F(x + 2p) — F(x — 2p) is an odd function with respect to p. □ 

To use the above symmetries we need the following lemma. 
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Lemma 1 If a smooth (analytic) function germ f : R, — > R, is not flat then the vector 
field germ f(u)d u on R, is equivalent to u k d u with respect to a certain smooth (analytic) 
coordinate transformation u = F(u), where k 6 N is such that d is the first non- 

vanishing derivative at . 



u=0 



Proof: F(u) must satisfy ODE f(u)F'(u) = (F(u)) k . The existence of F(u) with F'(u)\ u=0 ^ 
is easily verified in both smooth and analytic cases. □ 

Theorem 3 Suppose the solutions of h!3\) and those of 

a(x,y)dx + (3(x,y)dy = 0, (21) 

where a(x, y), f3(x, y) are non-flat functions at (0,0), form together a hexagonal 3-web. Then 
there is a local symmetry of [13]) at (0, 0) that maps equation [2l\) to one of the two following 
forms for y > 0: 

[(x + 2^y) k + {x- 2^y) k ]dx - -U(ar + 2^y) k - (x - 2^y) k ]dy = 0, or 

(22) 

y/y[(x + 2^y) k - (x - 2^y) k ]dx - [(x + 2^y) k + (x - 2^y) k ]dy = 0, 

where k is a non-negative integer. In particular, if a(x,y), /3(x,y) are non-flat functions with 
(a(Q, 0), (3(0, 0)) / (0,0), one gets three normal forms: 

a) dx = 0, b) dy = 0, c) 2dy - xdx = 0. (23) 

Moreover, if equation [2l\) for y > is equivalent to l23a) or 12%) then it can be reduced to 
I23\a). respectively $2 3b) by a symmetry of U3\) in some neighborhood of the point (0,0) 



Proof: Let us introduce operators U, V of differentiation along the curves of the foliations T and 
i(T) on M: 

U = 8 p + 2d x , V = 8 p -2d x . (24) 
Then these operators commute and satisfy the following relations: 

U(u) = 0, U{v) = 4, V(u) = -4, V{v) = 0. 

Consequently a direction field on M, whose integral curves form a hexagonal 3-web together 
with T and i(J~), must be generated by a vector field, commuting with U and V (for the details 
see [3], p. 17). Such a vector field has the form 

Y = f(v)U + g(u)V. 

The direction field generated by Y is the lift to M of the direction field induced by (|2ip iff 
Y A di(Y) = 0. This gives g = ±f. Projecting from M to the plane one obtains 

-U/(x + 2Jy) - f(x - 2^y)]d x + [f(x + 2jy) + f(x - 2^y)\d y 

for g = f and 

[f(x + 2^y) + f[x - 2^)]d x + y/y\f{x + 2jy) - f(x - 2^y)]d y , 
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for g = —f. Now applying symmetry (|17|) with F satisfying f(u)F'(u) = {F(u)) k (see Lemma 
[1]) we reduce f)21 1) to one of the forms (|22p for y > 0. If the found symmetry maps (|2ip to 
(123a ) than we can construct a diffeomorphism </? such that it is the identity for y > and maps 
integral curves of (|2ip to the lines x = const as follows. As equation (|2ip is not singular at 
(0, 0) it has a smooth first integral /(x, y) coinciding with x for y > 0. Then <^> is defined by 
(x,y) i— > (I(x,y),y). Similarly, for ([23b ) we define by (x,y) i— > (x,/(x,y), where I{x,y) is the 
first integral of (|2ip . coinciding with y for y > 0. 



2.2 The case of non-Legendrian criminant 

Theorem 4 Ze£ |7]) be an implicit ODE such that the corresponding surface M is smooth. 
Suppose its criminant C is a smooth curve, the projection ir : (x,y,p) — * (x,y) has a fold 
singularity at m £ C, and the contact plane is not tangent to M at m. Then {7J) is locally 
equivalent to 

p 2 = x (25) 

with respect to some diffeomorphism (p : U — > U, where U,U £l 2 are neighborhoods of (0,0), 
7t(to) and ^(0,0) = 7r(m). 

The proof is given in [I], p. 27. Similar to the case of Legendrian criminant, the diffeomorphism 
ip is not unique. 



Proposition 3 Equation (25\l has an infinite symmetry pseudogroup. Its transformations are 
given by 



x= \^-[F(3y + 2x^)-F(3y-2x^)} 2 , y = l{F{3y + 2^) + F(3y - 2xy/x)). (26) 
V 16 6 

Here F is a smooth function subjected to F'(u)\ u=0 ^ 0. Infinitesimal generators of this pseu- 
dogroup have the form 



{/(3y + 2xy^) - /(3y - 2x^)}8 x + {/(3y + 2x^) + /(3y - 2xy/x)}d y , (27) 



1 

\fx 

where f is an arbitrary smooth function. 

Proof: On the surface M defined by ([2]) we choose (y,p) as local coordinates. Solutions of (|25p 
define foliations T and by 

u := 3y — 2p 3 = const, v := 3y + 2p 3 = const, (28) 

where i is the involution i : M — ► M, (x, y,p) — > (x, y, —p). 

To prove the finite transformation formulas observe that the symmetry group transformation 
(x, y) I— s- (x, y) lifted to M must satisfy 

3y + 2p 3 = F{3y + 2p 3 ), 3y - 2p 3 = F(3y - 2p 3 ). 

In fact, to preserve the foliations T and it is necessary that 3y + 2p 3 = F(3y + 2p 3 ), 3y — 
2p 3 = G{3y - 2p 3 ). On the line p = one has F(3y) = G(3y) = 3y. This implies §2E§. The 
condition F'(u)\ u=0 ^ is equivalent to non-degeneracy of the Jacobian of the transformation. 
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Consider now an infinitesimal symmetry £(x, y)d x + w(x, y)d y . Lifted on M it turns to 

X = i](p 2 ,y)d y + g(y,p)d p . 

(On can write down an explicit expression for g, but we do not need it.) Since X is a symmetry 
of (|25p it must satisfy 

X{u) = 3/(u) = 3/(3?/ - 2p 3 ), X(t>) = 3k(y) = 3k(3y + 2p 3 ), 

for some smooth functions /, k. This is equivalent to 

v(p 2 ,y) = 2p 2 g{y,p) + f(3y - 2 P 3 ), v(p 2 , y) = -2p 2 g(y,p) + k{3y + 2p 3 ). 

Substituting p = into the difference of the above equations 

±P 2 g(y,p) = H3y + 2p 3 ) - f{3y - 2p 3 ) 

one gets k = f. Hence the functions r\ and g are well defined by 

V(P 2 ,V) = \(f{Zy + 2p 3 ) + f(3y - 2p 3 )), g(y,p) = ^(/(3y + 2p 3 ) - f(3y - 2p 3 )). 

Note that g(y, 0) = 0, hence X is tangent to the criminant and therefore lowerable (see [2]). Up 
to scaling the lowered operator X becomes (|27p . The extension of the defined transformation 
for x < is again justified by Malgrange's Preparation Theorem. (See the detail in the proof of 
Theorem [2j) 

Theorem 5 Suppose the solutions of {25\) and those of 

a(x,y)dx + (3(x,y)dy, (29) 

where a(x, y), (3(x, y) are non-flat functions at (0,0), form together a hexagonal 3-web. Then 
there is a local symmetry of \25\) that maps equation \29\) to one of the two following forms for 
x > 0: 

]=[(3y + 2x^) k - (3y - 2xy/x) h ]dx - \fl6[{3y + 2x^x) k + (3y - 2xy/x) k ]dy = or 



\{3y + 2x^/x) k + (3y - 2xy/x) k ]dx - \^16^/x[(3y + 2xJ£) k - (3y - 2xy/x) k ]dy = 0, 

where k is non-negative integer. In particular, if a(x,y), f3(x,y) are non-flat functions with 
(a(0, 0), /3(0, 0)) 7^ (0,0), then equation \29\) can be reduced to one of the following two normal 
forms in some neighborhood of the point (0,0): 

a) dx = 0, b) dy = 0. (30) 
Proof: Let us introduce operators U, V of differentiation along the curves of the foliations T and 

U = d p + 2p 2 d y , V = d p -2p 2 d y . (31) 
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Then the operators and commute and satisfy the following relations: 
\u(u) = 0, \u(v) = 12, \v(u) = -12, \v{v) = 0. 

'P P p 

A direction field on M, whose integral curves form a hexagonal 3- web together with T and i(J-), 
must be generated by the vector field, commuting with -^U and ^zV. Such a vector field has 
the form 

Y = f(v)\u + g{u)\v. 

The direction field generated by Y is the lift to M of the direction field induced by (|29p iff 
Y A di(Y) = 0. This gives g = db/ (compare with the proof of Theorem [3]). Projecting from M 
to (x, y)-plane one obtains 

[f(3y + 2x^) + f(3y - 2x^)]8 x + y/E[v@y + 2x^) - /(3y - 2x^)]d y 

for g = f and 

i[/(3y + 2x^) - /(3y - 2xy/x)]d x + [v(3y + 2xjx) + /(3y - 2x^)]a y 



for for g = —f ■ Now applying symmetry (|26l) with i 7 satisfying f(u)F'(u) = (F(u)) k (see 
Lemma [1]) we complete the proof. The details can be found in the proof of Theorem [3j 

Remark 1 . Real analytic versions of Theorems [3] and [5] are true without the stipulations y > 
and x > respectively. 

Remark 2. One can not extend the claim of Theorem [3] for y < in the smooth case for 
equation (|23b ). Its solutions are parabolas y = x 2 /4 + C . They cross the line y = in two 
points if C < 0. If one smoothly deforms equation (|23|) in the domain y < then the solution 
of the deformed equation starting from some point (xq, 0) with xq < will not necessarily pass 
through (—xo,0), i.e. this solution returns to a "wrong parabola". 

Corollary 1 If the following conditions hold for implicit ODE at a point m = (xq, yo,Po) £ 
M: 

1 ) ODE (OJj has a hexagonal 3-web of solutions, 

2) po is a double root of HP at ir(m) = (xo,yo), 

3) regularity condition is satisfied, i.e. rank((x,y,p) i— ► (F,F p ))\ m = 2, 

then its criminant is either Legendrian or transverse to the contact plane field in some neigh- 
borhood of m 

Proof: Denote by Ct the closed set of points on the criminant C, where the contact plane is 
tangent to C . Suppose m is not a point of C \ Ct and not an interior point of Ct- Then m is 
a boundary point of Ct- Now Theorem [5] implies that for each point m! sufficiently close to m 
and such that m! ^ m, m! ^ Ct equation (TjQ) is locally equivalent to a product of an explicit 
ODE dx = and quadratic equation (|25j) . i.e. the solutions of the linear factor are tangent to 
the discriminant curve at ir(m') and therefore at 7r(m). Further Theorem [3] implies that for each 
point m' sufficiently close to m and such that ml ^ m, ml £ Ct equation (TIJ is locally equiva- 
lent to a product of an explicit ODE dx = and quadratic equation ([TBI) , i.e. the solutions of 
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the quadratic factor are tangent to the discriminant curve at 7r(m') and therefore at 7r(m). But 
that means that the root po is triple. Thus our assumption is false and the corollary is proved. □ 

Remark 3. As follows from the above proof the hypothesis of Corollary [1] also implies that 
there is no isolated points of tangency of the contact plane and the criminant. 



3 Normal form for an ordinary cusp point 

In this section we use the results of the previous one to establish normal forms for the case of a 
cusp singularity of the projection ir on M. Regularity condition for a triple root pq implies 
immediately that the projection ir has a cusp point at m = (xo,yo,Po) £ M. We start with a 
Legendrian criminant, then consider non-Legendrian criminant and finally show that one can 
not "glue" Legendrian criminant with non-Legendrian one at the cusp point. 

Lemma 2 If the following conditions hold for implicit ODE (Op at a point m = {xq, yo,po) £ M: 

1 ) ODE (C]j has a hexagonal 3-web of solutions, 

2) pq is the triple root of |7P at ir(m) = (xo,yo)> 

3) regularity condition is satisfied, i.e. rank((x,y,p) i— ► (F,F p ))\ m = 2, 
then it is locally equivalent to 

p 3 + A{x,y)p + B{x,y) = 0, (32) 

where 

1) the projection ir has an ordinary cusp singularity at (0,0,0) with A(0,0) = B(0,0) = 0, 



d(A,B) 

i isvuiuuouoco Utb yu,vj, i.e. 

3) B x {0,0) =0. 



2) A, B are local coordinates at (0,0), i.e. d ^ x y ^ 



7^0 

x=y=0 



Proof: Since equation (pQ) has the triple root po at n(m) it is locally equivalent to some cubic 
equation Q. Further, the coefficient by p 2 in this cubic equation is killed by a coordinate 
transform of the form y = f(x,y), x = x, satisfying 

3fx{x,y) + a(x,y) = 0. 

This transform respects the regularity conditions. Thus our implicit equation F = becomes 

F(x, y,p) = p 3 + A(x, y)p + B(x, y) = 0. 

Without loss of generality it can be assumed that (xo,yo) = (0,0). As the equation p^ + 
A(0, 0)po + B(0, 0) = has a triple root holds po = 0. Therefore the functions A, B must also 
vanish at (0,0). Now regularity condition J4} at m = (0,0,0) 6 C reads as 



| A xP + B x A y p + By 3p 2 + A 

A x A y 6p 



B x By 

\ A-r A.*/ 
x = y=p=0 \ x y 



= 2. 

x=y=p=0 



Thus claims 1) and 2) are proved. Moreover, the discriminant curve A = ir({(x,y,p) : p 3 + 
A(x, y)p + B(x, y) = 3p 2 + A(x, y) = 0}) has an ordinary cusp at 7r(m) = (xq, yo). 
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If solutions of equation (|32|) form a hexagonal 3-web the curvature of this 3-web must vanish 
identically. This is equivalent to the following cumbersome partial differential equation for the 
functions A, B: 

(4A 3 + 27B 2 )(9BA XX - 2A 2 A xy + 6ABA yy - 6AB XX - 9BB xy - AA 2 B yy )+ 
+W8A 2 BA x B y - W8AB 2 A x A y + \Q2B 3 A 2 y + 40A 4 A y B y - 108A 2 BA 2 X + 
+216A 2 BB 2 - 36A 3 B x B y + 108A 2 BA y B x - 378AB 2 A y B y - 405B 2 A X B X + ^ ' 
-48A 3 BAl + 8A 4 A x A y + 2A3B 2 B x B y + 8AA 3 A X B X + 32AABB 2 = 0. 

This equation is obtained by direct lengthy but straightforward computation. (Expressions for 
the corresponding web curvature for a cubic ODE can be also found in [19] and |22j). As was 
shown above the functions A, B can be taken as local coordinates around (0, 0). Then all partial 
derivatives of A, B with respect to x and y are smooth functions of A, B. The homogeneous 
part of second order of Taylor expansion of l.h.s. of (|33|) around (0, 0) is 

-405B 2 A X (0, 0)B X (0, 0) + 243B 2 B X (0, 0)B y (0, 0) + 324AB(B X (0, 0)) 2 . 

It must vanish. In particular, B x (0, 0) 2 = as the coefficient by AB. 



3.1 The case of Legendrian criminant 

Theorem 6 If the following conditions hold for implicit ODE at a point m = (xo,yo,po) £ 
M: 

1 ) ODE |7]j has a hexagonal 3-web of solutions, 

2) po is the triple root of |2P at ir(m) = (xQ,yo), 

3) its criminant C is a Legendrian curve with rank((x, y,p) *— > (F,F p ))\c = 2, 
then it is locally equivalent to the following Clairaut equation 

P 3 + PX -Y = 0. (34) 

Proof: Lemma [2] reduces equation © to (EH) with A(0,0) = £(0,0) = ^(0,0) = and 
d af' B ) 7^ 0. Thus the tangent plane T m M to the surface M at m = (0, 0, 0) is the 

o\x,y) x=y=0 

plane y = and the tangent line to the discriminant curve at (0, 0) is y = 0. (The condition 
B x (0, 0) = for the case of Legendrian criminant can be obtain also by the following geometrical 
consideration. By Proposition[T]the characteristic field r, given by ([5]), can be smoothly extended 
to the criminant C. As r is transverse to C the projection of integral curves of r are smooth 
curves in M 2 tangent to the discriminant curve A at the origin (0, 0). As po = this implies the 
claim.) With B x (0,0) = one gets from the regularity condition 

4,(0,0)^0, ^(0,0)^0. 

Therefore one can choose p, A as local coordinates on the surface Matm= (0, 0, 0) and A, B 
as local coordinates on M 2 ,0. In these coordinates the projection it is the Whitney map. The 
criminant is parameterized by p as follows 

A = -3p 2 , B = 2p 3 . 
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Its projection is the discriminant curve A := {(A, B) : 27B 2 + AA 3 = 0}. The set of points 
projected to the discriminant curve is the criminant itself and the following curve 

D:={(p,A) GM: A = ~p 2 }. (35) 

This follows from the observation that the value — 2p is the third root of (|32p at the discriminant, 
where the double root is p. The curve D is tangent to C at 0, thus the characteristic field r is 
transverse also to D. 

Consider the following map 

/:(R 2 ,0)^(M,0), (p,q)^(p,A) = (p,-q 2 -pq-p 2 ). (36) 

This map has a fold singularity on the line L\ := {(p, q) : p + 2q = 0}. This line is mapped 
by / to the curve D since — q 2 — pq — p 2 = — [q + |) 2 — \p 2 ■ Note that if f(p, q) = (p, A) then 
f^ 1 (p, A) = {(p, q) U (p, — p — q)}. The pull back f of the characteristic field r by /* from M, 
to R 2 , must be tangent to the kernel of df ', i.e. to the vector field <9q, since r is transverse to 
D. Moreover, the foliation T\ of integral curves of f is invariant with respect to the following 
linear involution 

51 :R 2 ^R 2 , (p,q)^(p,-p-q). 
Consider also the following two linear involutions 

5 2 :R 2 ^R 2 , (p,q) » (-p- q,q), g 3 : R 2 ^ R 2 , (p,q) ^ (q,p). (37) 

The linear maps 51,52,53 generate the group B3, the symmetry group of equilateral triangle, 
which can be viewed as the group of linear transformations of the plane p + g + r = 0inR 3 , 
generated by permutations of the coordinates (p,q,r) in R 3 . The orbit of a point (p,q) under 
this group action is the inverse image of the point (A,B) = n(f(p,q)) € R 2 under the Vieta 
map V := 7T o /. Therefore the three foliations !F\, JT 2 := g% (^1) and T3 := 53(^1) form a 
hexagonal 3-web. Moreover, this 3-web is not singular at (0, 0) and has the symmetry group ID3 
generated by {51,32,53}- Note that for Clairaut equation ([M]) the above defined three foliations 
are p = const, p + q = const and q = const respectively. 

Now we are ready to construct the diffeomorphism ip that transforms the given ODE to 
normal form (|34|) . Consider a domain U, C R 2 , such that the 3-web formed by the foliations 
JFi, T2 and Tj, is regular in U. Let 71 be the integral curve of f that passes through the origin. 
Pick up a point u = (pi, q±) G 71 on this curve and draw the Briangon hexagon around through 
u. (Let us recall the construction of the Briangon hexagon: one draws three curves j{, i = 1, 2, 3 
of the foliations i = 1, 2, 3 through the origin, picks up a point on one of this curves, say 71, 
and then goes around the origin along the foliation curves, swapping the family whenever one 
meets one of the 7$. The web is hexagonal iff one gets a closed hexagonal figure for any choice of 
the central "origin" point and u. See Fig. [3] on the left.) Let us choose u so that the following 
conditions hold: 

1) 51 > 0, 

2) the Briangon hexagon around through u is contained in U. 

Then there is a unique local homeomorphism ip : R 2 , — > R 2 , such that 

1) V(0,1) =u, 

2) VU,0) =g 3 (u), 
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Figure 3: Briancon's hexagons and their inverse images under the diffeomorphism ip with u = 

if;(v). 



3) it maps the foliations p = const, p + q = const and q = const to the foliations J-\, T% and T-$ 
respectively. (See Fig. [3|). In fact, the points u and g^ (u) lies on the same curve of the foliation 
J-2 since the involution g^ is a symmetry of Further, there is a unique diffeomorphism, 
mapping the triangle (0, 0), (0, 1), (1, 0) to the "triangle" (0,0), u,g^{u) formed by the curves 
of the foliations Ti and JF 3 (see [3] p. 15). This map is uniquely extended to the whole 
hexagon. Moreover, the constructed homeomorphism ip is equivariant with respect to the action 
of B3 defined above. The map ip is smooth (analytic) if the foliations T\,TiiJ~z are smooth 
(analytic). Really, according to [3] p. 155, there exists a smooth (analytic) map, taking the 
foliations p = const, p + q = const and q = const to T\, T2 and respectively, and this map is 
uniquely defined by specifying the inverse image of u. Thus this map should coincide with the 
above homeomorphism ip. 
Now consider the map: 

7T o / o V : M 2 , - R 2 , 0, (P, Q) 1 * (A, B) = tt(/(^(P, Q))). 

In coordinates it reads as 

A = a{P,Q), B = (3(P,Q). 

Observe that the above map is symmetric with respect to the action of B3. Then by the results 
on smooth functions, invariant with respect to finite group action, the functions a and (3 must 
depend only on the basic invariants of the above B3 group action (see [E]) and [23]): 

A = a(X,Y), B = (3(X,Y), 

where 

X = -P 2 -PQ-Q 2 , Y = PQ{P + Q). 

We claim that a, f3 are the components of the searched for diffeomorphism </? 
consider the following commutative diagram: 

R 2 )0 R 2 )0 

I V [n of (39) 

M 2 ,0 R 2 ,0 



(38) 

To prove that 
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where V is Vieta's map (|38p . Applying the same results on symmetric functions to 8 := V oip' 1 
we see that the differentiable map 6 is inverse to cp inside the "cusped" domain where our ODE 
has 3 distinct real solutions. This completes the proof in the real analytic case. For the smooth 
case we apply tp and for the reduced equation we consider the first integral of the direction field 
r that coincides with p on the part M3 of M that is projected to the domain with three real 
roots. Further, it is easy to construct through homotopy the 7r-lowerable diffeomorphism <p' of 
M that is identity on M3 and moves the integral curves of r to that of (134j) . The searched for 
diffeomorphism is ip' o ip 



3.2 The case of non-Legendrian criminant 

Theorem 7 // the following conditions hold for implicit ODE {T]j at a point m = (xq,pq,pq) £ 
M: 

1 ) ODE (OJj has a hexagonal 3-web of solutions, 

2) po is the triple root of |7]) at 7r(m) = (xo,yo), 

3) the criminant C is transverse to the contact plane field in some punctured neighborhood of m 
and rank((x, y,p) i-> (F, F p ))\ c = 2, 

then it is locally equivalent to 

P 3 + 2PX + Y = 0, (40) 

within the domain, where has three real roots, if F is smooth, 
and in some neighborhood of (0,0), if F is real analytic. 

Proof: We follow the proof scheme for Theorem [6j Namely we consider the pull-back of the form 
dy—pdx to M 2 , by the Vieta map vro/, where / is defined by (f36l) . duplicate this pull-back form 
by linear involutions ()37[) and find a local diffeomorphism of M? , matching " lifted" 3-web of our 
equation and that of (|40p . The difference to the previous case of Legendrian criminant is that 
now the web is singular; each foliation Ti has a saddle singular point at f~ l (m). Therefore the 
classical results on hexagonal 3-web are not of much use to find the diffeomorphism " upstairs" . 
We construct it through a homotopy of the first integrals of the corresponding foliations. 
•Differential forms of the foliations. 

By Lemma [2] equation {TJ is equivalent to ([32]) with A(0, 0) = B{0, 0) = B x (0, 0) = 0, §^ / 
and 

4,(0,0)^0, ^(0,0)^0. 

Thus the tangent plane T m M to the surface M at m = (0, 0, 0) is the plane y = and the tangent 
line to the discriminant curve at (0, 0) is y = 0. Therefore one can choose p, A as local coordinates 
on the surface M, in and A, B as local coordinates on R 2 , 0. Thus x = X(A, B), y = Y(A, B), 
where A = —p 2 — q 2 — pq and p,q,r = —p — q are roots of (|32j) . As easily follows from Theorem 
[5j the kernel of the pull-back form tt* (dy — pdx) is tangent to the curve D defined by (|35p. That 
means that the kernel of the form 

f*(7T*(dy-pdx)) = (2p + q)(-Y A +pX A + qY B -pqX B )dp+(2q+p)(-Y A +p(X A + Y B )-p 2 X b )dq 

is tangent to the line L\ := {(p,q) : p + 2q = 0}. Writing the above form as uj\ := (2p + 
q)P(j), q)dp + (2q + p)Q(p, q)dq with suitable P, Q and passing to the coordinates q, s = 2q + p 
one gets 

uji = (2s - 3q)Pds + (sQ - 2(2s - 3q)P)dq 
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hence the tangency condition implies P\ s =o = 0. By the Hadamard lemma P = sP hence one 
obtains 

u x = s{(2s - 3q)Pds + (Q - 2(2s - 3q)P)dq}. 

Now consider the expression for P = —Ya + (s — 2q)XA + qYs — (s — 2q)qXs- One has Ya = 0, 
Xa 7^ from A x B y ^ 0, B x = 0. Since A is quadratic and B is cubic in p, q, the term Ya does 
not have linear terms in p, q. This implies -P(0, 0) ^ 0. Using again the condition on L\ one 
obtains 

Q - 2(2s - 3q)P = sQ. 

Let as normalize the forms vanishing each on its own family of solutions to satisfy a\ + 02 + 
(73 = 0: 

o"i = (g — r){dy — pdx), 02 = (r — p){dy — qdx), 03 = (p — q)(dy — rdx). (41) 

As shown above the pull-back of u\ is 

(fi = (2g + p) 2 {(2p + g)Pdp+(2(2p + g)P + (2g+p)Q)dg} (42) 

• Connection form. 

Following [4j consider the area form 

n := ox A o 2 = 02 A (fa = (f 3 A ai = (Fa^b - Y B X A ){p - q) 2 (2p + ^) 2 (2g + pfdp A 

and the connection form 

7 := h 2 d 1 - hid 2 = h 3 (J2 - h 2 <T 3 = hid 3 - h 3 di, 

where hi are defined by 

ddi = hi^l. 

Using d32J) on obtains by direct calculation that 

, ^1 

(2p + q) 2 (p - q) 2 ' 

where R\ is a smooth function of p, q. Applying the cyclic permutation p — * q, q — > r, r — > p, 
one gets: d 2 = (2p + q) 2 d 2 , h 2 = ^ 2 q+pf^(p-q)' :! > wnere ^2 an d d 2 are smooth. Therefore 

= 1 
(p - q) 2 ' 

with a smooth form 7. Observe that the connection form 7 is symmetric with respect to the 
linear transformation group O3, generated by gi,g 2 ,93- Therefore 7 is smooth, i.e. (p — q) 2 
divides 7. 

•Existence of first integrals. 

As the web is hexagonal the connection form 7 is closed. Therefore there exists a unique B3- 
symmetric function /j,, satisfying 

d/i = //7, p(0,0) = l. 
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Further, the forms fidi are also closed, thus defying functions U( by 

d(ui) = fidi, Ui(0,0) = 0, 

satisfying the following equation, which is equivalent to the hexagonality of the web: 

m + u 2 + u 3 = 0. (43) 

Observe that the function u\ is skew-symmetric with respect to g\\ 

g* 1 (u 1 ) = -u 1 . (44) 

This follows from (|4ip and from the invariance of /i. Applying Hadamard's trick one estimates 
u\ as follows: 



ui(p,q) 



—ui(tp,tq)dt = J f p— U!(tp,tq) + q—ui{tp, tq) \ dt = 

= (2q + pf J n(tp, tq) (p(2p + q)P(tp, tq) + q(2(2p + q)P(tp, tq) + (2q + p)Q(tp, tq)fj t 3 dt. 
Collecting similar terms, using P(0, 0) ^ and integrating one has 

«l(p,g) = (2q + p) 3 ((2p + q)P(p,q) +qQ(p,q)), where P(0, 0) / 0. (45) 
• Properties of the first integrals. 

It follows from Malgrange's Preparation Theorem that any smooth function F of (p, q) can be 
represented in the form 

F(p,q) = F (A,B) + P F 1 (A,B)+qF 2 (A,B) + pqF 3 {A, B) + q 2 F 4 (A, B) +pq 2 F 5 (A,B), (46) 

with smooth functions i^. In fact, the identities p 2 = —pq — q 2 — A, p 3 = — pA — B, p 2 q = 
-pq 2 + B, q 3 = -qA-B imply (p, q) 4 C (A,B) and £(R 2 )/(A, B) = R{l,p,q,pq,q 2 ,pq 2 }. (Here 
<S(M 2 ) is the local algebra of smooth map germs at (0,0), (p,q) its maximal ideal, generated 
by the coordinate functions p and q, A : (p,q) i— ► — p 2 — pq — q 2 , B : (p,q) i— ► p 2 q + q 2 p, 
and M{l,p, q,pq, q 2 ,pq 2 } is the real vector subspace of £(M. 2 ), spanned by l,p,q,pq,q 2 ,pq 2 -) 
Moreover, inside the "cusped" domain with 3 real distinct solutions of our ODE the functions 
Fi are uniquely determined by F. For F = u\ property (|44p implies F^ = — F3, F 2 = 2F\ — AF§, 
F = — AF3 — BF 5 /2. Applying — g% and —g^ to u\ one gets the other two first integrals u 2 
and U3, whose representations in form (|46p are easily read from the representation of u\. Now 
identity (g5D implies F 5 = 0. Thus 

«i(p,g) = (2g + p)(Fi(A,5) +pF 3 (^,B)). (47) 

Using (|43|) one can write 

Fi(AS)+pF 3 (AS) = (2g + p) 2 G(p,g). 
Representing the function G as 

G(p, g) = G (A, B) + pdiA, B) + qG 2 (A, B) + pqG 3 (A, B) + q 2 G i (A, B) + pq 2 G 5 (A, B) 
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and substituting this representation into the above equation, one obtains 

G(p, q) = ^AG A (A, B) + pG^A, B) + pqG A (A, B) + q 2 G A (A, B). (48) 
•Equivariant homotopy. 

For equation (|4"0j) the first integral u\ is uq(j>, q) := p(2q+p) 3 . Let us scale u\ so that P(0, 0) = \ 
in (|45p . We claim that the family of functions 

itt(p,g) := u (p,q) +t(ui(p,q) -u (p,q)) 

is equivariantly 7£-trivial, i.e. for any t £ [0, 1] there is a diffeomorphism tpt, equivariant with 
respect to above B3 group action, such that 

u t oip t = u . 

To prove this it is enough to find B3-equivariant vector field (,(p,q,t)d p + r](p,q,t)d g satisfying 
the following homotopy equation: 

dut(p,q) , du t (p,q) du t {p,q) 
£(p, q, t) — + i](p, q, t) — + — = 0, f (0, 0, t) = 77(0, 0, t) = 0. 

A general form of a D3-equivariant vector field is given by 

C=pa(A,B)+(^+pq + q 2 )(3(A,B), r\ = qa(A, B) — (\a + q 2 \ (3(A, B), (49) 



(see [2] or derive it from the representations of £,77 in form (|46p ). Observe that the difference 
(«i(p,g) - uo(p,q)) = also has form (08]): 

ui(p, q) - uo(p, q) = ^AL(A, B) + P K(A, B) + pqL(A, B) + q 2 L(A, B) 

with K(0, 0) = due to the chosen scaling of u\. Solving the homotopy equation yields the 
following expressions for a and /?: 

6K + (6K 2 + 2A 2 (K B L - KL B ) + 9B(KL A - K A L) + WAL 2 )t 
a = 



M 

12L + (6A(KL A - K A L) + 9B(KL B - K B L) + 3KL)t 



M 
where 

M = —24 + [-48K - Y1AK A - 18B(K B + 2L A ) + 8A 2 L B ]t+ 

+ [-24K 2 - 2A(6KK A + 25L 2 ) + 3B(15K A L - 6KK B - 12KL A )+ 

+2A 2 {AKL B - WLL A - 5K B L) - 30ABLL B + (27B 2 + 4A 3 )(K A L B - K B L A )]t 2 . 

M does not vanish at (0, 0) since K(0, 0) = 0. The claim on 7£-triviality of the family ut(p, q) is 
proved. 

•Diffeomorphism. 

We have proved that the diffeomorphism 

ip := 
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maps the fibres of uq, i.e. the curves {(p,q) : u (p,q) = const} to those of u\. Therefore, being 
equivariant, (p maps the foliations T% of (|40p to that of our equation (|32p. Now diagram (|39|) 
defines again the desired diffeomorpfism ip. □ 



Remark. The pictures of the B3-symmetric hexagonal 3-web, defined by the solutions of (|40p 
and lifted to the plane p + q + r = 0, is presented in Fig. [J] on the left. It consists of 3 foliations, 
one of them is shown in Fig. |4] in the center. On the right is the fundamental domain of 
ID>3-group (compare with Fig. Q] on the right). The flower- like form on the left suggests that 
the web is actually symmetric with respect to the symmetry group ©6 of regular hexagon. In 
fact, it is the case since the fibers of the first integrals are permuted by the following symmetry 
g 4 :R 2 ^R 2 , (p, q) h- (-p, -q). 




Figure 4: D3-symmetric hexagonal web of 3 foliations with saddle singularities (left), one of the 
foliations (center) and the fundamental domain of D3 (right). 



4 Proof of the classification theorem 

Now we can prove Theorem [TJ If a point m = (xo,yo,po) G M is regular then our equation is 
locally equivalent to (fl~2"V ) by definition. If po is a double root then the regularity condition Q 
implies that m is a fold point of the projection ir. Futher Corollary [1] implies that the criminant 
is either Legendrian or transverse to the contact plane field in some neighborhood of m. Thus 
by Theorems [3] and [5] the equation is locally equivalent either to (|12t ii) or to (|12t v). Finally if po 
is a triple root and the criminant is either Legendrian or transverse to the contact plane field in 
some punctured neighborhood of m then by Theorems [U] and [7] the equation is locally equivalent 
either to (|12t ) or to (|12t i). To complete the proof we show that Legendrian and non-Legendrian 
parts of the criminant can not be glued together at a cusp point. By Lemma [2] our equation is 
equivalent to ([32]) with (xo;2/0;Po) = (0,0,0), B x (0,0) = 0. Suppose the criminant is transverse 
to the contact plane field for p > and Legemdrian for p < 0. For any point m! 7^ m with 
p > on the curve D defined by (|35p the direction field r is tangent to D by theorem This 
condition reads as 



3 

^(^4 + -p 2 ) A {dy - pdx) 



= 0. 

D 
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In coordinates p, x on M it can be rewritten as follows 

= 0. 



3 

(dA + -pdp) A ((3p 2 + A)dp + [B x + {A x + B y )p + A y p 2 )dx) 



D 



Substituting A = — \p 2 and 



dA = A x dx + Aydy = A x dx - - v ((3p 2 + A)dp + (B x + A x p)dx) 

By + AyP 

into this equation one gets 
3 

-p 2 (A X By + Ay(A X + By)p + - P ( B a + (4, + By)? + Aytf){By + 



= 0. 

D 



Parameterizing the curve Z? by p, expanding the above equation by Tailor formula at p = and 
equating the coefficient by p 2 to one obtains 



3 

By(-A X -(A X +By)) 



= 0, 

x=y=0 



which implies 

253,(0,0) -4,(0,0) =0, 

since B y (0,0) / 0. (We have used the Taylor formula B x = B xA A+B xB B + ... = B xA (-3/4p 2 ) + 
B xB (-l/ 4p 3 ) + ....) 

On the other hand, for any point mf ^ m with p < the contact form vanishes on the 
criminant: 

dy — pdx\ c = 0. 

In coordinates p,x on M it can be rewritten as follows 

(B x + (A x + B y )p + A yP 2 )\ c = 0. 

Now the Tailor expansion at p = for C parameterized by p gives A x + B y = (B x does not 
have linear in p terms). Comparing with the condition above on the non-Legendrian part one 
gets Ac (0,0) = and therefore 5^(0,0) = which contradicts Lemma [2l ^ 

Remark. Unfortunately, the annoying stipulation in Theorem []] for the smooth case ()12t ) 
can not be omitted to guarantee the existence of the diffeomorphism ip reducing ODE under 
consideration to (|12t ) in some neighborhood of m if one stays within the framework of geometric 
Definition [TJ A necessary condition for that is the existence of the first integral of r in the form 
f{p-,x) 2 g{p,q)' i = (p* (p 2 (x + 3p 2 / '8) 3 ) . Here (p is the lift to M of the searched for diffeomorphism 
and p 2 (x + 3p 2 /8) 3 is the first integral of r for (1120 . It is not hard to find a counterexample 
which does not have such an integral in the form f{p,x) 2 g(p,q) 3 with not vanishing df,dg at 
(0,0). This drawback is repaired as follows. One replace definition [T] with a less geometric one. 

Definition 2 We say that implicit ODE \3ty has a hexagonal 3-web of solutions if A, B satisfy 
PDE \33\) and the domain, where has 3 real roots Pi,P2,P3 is not empty. 

The proof of Theorem [7] is easily modified for the case of one real root pi = p and two complex 
conjugated roots p2,3 = — p/2 ± iz. The form a\ turns out to be pure imaginary but the 
connection form 7 is real. All analytical properties being the same, one finds the diffeomorpfism 
<p similarly through homotopy. 
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5 Concluding remarks 



• Symmetries of the normal forms. The solutions of equation (|12| y) are the lines dx = 0, 
dy = and dx + dy = 0. The symmetry group of (TT2V ) is generated by the following operators: 

Xi = c\., X 2 = <9 a , X 3 = x5 x . + yd y . 

Thus the symmetry pseudogroup of a cubic implicit ODE with hexagonal 3- web of solutions is at 
most 3-dimensional. In a neighborhood of the projection of a regular point m G M it is generated 
by the above three operators Xj in suitable coordinates. The coordinate change becomes singular 
on the discriminant curve and not all symmetry operators "survive" at 7r(m) £ A. The symmetry 
pseudogroups of equations (fT2l ii) and (fT2l v) at a fold point are generated by 

Yx = xd x + 2yd y , Y 2 = d x 

and 

Yx = 2xd x + 3yd y , Y 2 = d y , 

respectively. This easily follows from Propositions [2] and [3j Irreducible equations (|120 and p21i) 
have only one-dimensional symmetry pseudogroup at (0,0): 

Z = 2xd x + 3yd y . 

• Analytic properties. All equations in the given normal forms are integrable in elementary 
functions. 

• Implicit cubic ODEs with singular surfaces M. Suppose our cubic ODE factors out to 
3 linear in p terms p — fi(x,y) such that 2 of 3 smooth surfaces Mj := {(x,y,p) : p = fi(x,y)} 
intersect transversally along a non-singular curve, the solutions of these 2 factors being transverse 
to the curve projection into the plane. Then one can bring these two factors to the forms p = 
and p = 2x respectively. The symmetry pseudogroup of the quadratic ODE p(p — 2x) = is 
y = F(y), x = y/F{y) — F(y — x 2 ). If our cubic equation has a hexagonal 3- web of solutions 
then its third factor is generated by the vector field (a(y — x 2 ) + j3(y))dx + 2xj3(y)d y . As the 
functions a,/3 are arbitrary we can hope to "kill" only one of them by the above mentioned 
symmetry. Thus a general classification of all cubic ODEs will have functional moduli even if 
one impose hexagonality condition. (Note that if the third family of solutions in the example is 
transverse to the first two, we have (3 = —a and one gets a finite classification list.) 

• Other examples. The proof of Theorem [7] suggests the following procedure to generate cubic 
ODEs with a hexagonal 3- web of solutions: start with a function F(p, q) written in form (|46p 
with A = p 2 +pq + q 2 , B = pq(p + q), define G(p, q) := F(q, —p — q), H(p,q) := F(—p - q,p) 
and solve the following equations for Ff. 

gl(F)=±F, F + G + H = 0. (50) 

This gives four of six coefficients F{ as linear combinations of the remaining two " free" functions 
of A, B. Then the fibers of F, G, H define a hexagonal 3-web, symmetric with respect to D3- 
group action, generated by the involutions gi,g2,93- The image of this web under Vieta map 
(p,q) l— * (p 2 + PQ + Q 2 > PQ(p + l)) is a hexagonal 3-web of solutions of some implicit cubic ODE. 
For example, starting with F = p — q one gets the following equation: 

2 1 

VP 3 - + xyp + — (2x 3 - 27y 2 ) = 0. 
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The solution 3-web of this equation is dual to that of (|12t ) (see [22J). Its surface M is not smooth 
at (0,0,0). Note that one can also start with F such that the "free" coefficients Fi have poles 
at 0. This approach linearizes the problem of finding local solutions of nonlinear PDE (I33|) . On 
the space of functions F satisfying ([50]) acts the pseudogroup of B3-equivariant transformations 
with the tangent space generated by vector fields defined by (|49|) . General classification of such 
functions with this equivalence group seems rather unpromising since the orbit codimension 
quickly becomes infinite. 
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